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A numerical simulation of silica aerogels is performed using diffusion-limited cluster-cluster ag- 
gregation of spheres inside a cubic box (with periodic boundary conditions) . The volume fraction 
c is taken to be sufficiently large to get a gel structure at the end of the process. In the case 
of monodisperse spheres, the wavevector dependent scattered intensity I{q) is calculated from the 
product of the form factor P{q) of a sphere by the structure factor S{q), which is related to the 
Fourier transform of g{r) — 1, where g{r) is the pair correlation function between sphere centers. 
The structure factor S{q) exhibits large-g damped oscillations characteristics of the short range 
(intra-aggregate) correlations between spheres. These oscillations influence the I{q) curve in the 
g-region between the fractal regime and the Porod regime and quantitative comparisons are made 
with experiments on colloidal aerogels. Moreover, at small-g values, S{q) goes through a maximum 
characteristic of large range (inter-aggregate) correlations. Quantitative fits of the maximum in 
the experimental I{q) curves of base-catalyzed aerogels are presented. In the case of polydisperse 
spheres, I{q) is calculated directly from a single aggregate simulation. It is shown that increasing 
polydispersity shifts the location of the cross-over between the fractal and Porod regimes towards 
low g-value. 



1. Introduction 

As revealed by Small Angle X-rays Scattering (SAXS) 
or Small Angle Neutron (SANS) experiments, Silica 
Aerogels [0-|3j are made of a disordered, but homoge- 
neous, array of connected fractal clusters which result 
from the aggregation of primary particles Q. In the 
case of colloidal aerogel the form factor of the individual 
particles is well defined and one can extract the precise 
form of the structure factor S{q) from the scattered in- 
tensity I{q). The analysis of the wave vector dependence 
of S{q), and I{q), has permitted the determination of two 
characteristic length scales which are the average size, a, 
of the particles and the average size, ^, of the clusters. 
Three distinct domains of wavevectors can be identified 
in both functions. At large q (q > a~^) S{q), and I{q), 
exhibits large-g damped oscillations, and the Porod law 
{I{q) ^ 9 respectively. At intermediate values of q, 
for < q < a~^, the fractal nature of intra-cluster 
particle correlations is revealed by a power law behavior 
(7^^, where D is the fractal dimension |^ of the clusters. 



At last, at small q values, for q « ^ , the scattering 
function saturate and eventually decreases as q tends to 
zero. 

Earlier studies ||^-[l0[| of the damped oscillations in the 
S{q) curve of colloidal aerogels have not been analyzed 
in terms of their all short range (intra-aggregate) fea- 
tures. On the other hand, the presence of a maximum 
in some S{q) and I{q) curves was puzzling for some au- 
thors who tried to fit the experimental data by 
considering single aggregate theories like the formula due 
to Fisher and Burford in which ^ plays the role of 
a radius of gyration. Others P,pp||T^ have used a semi- 
empirical formula obtained by introducing a cut-off func- 
tion to limit the fractal scaling, in which ^ enters as a 



correlation length. In all cases, these approaches do not 
take into account the inter-cluster interactions which give 
rise to a maximum in the scattering function. As far as 
long range properties are concerned, the single-aggregate 
approach is only valid for extremely diluted solutions of 
aggregated particles where the mean inter-aggregate dis- 
tance £ is much larger than the mean radius of gyration R 
of the aggregates. Otherwise, the theoretical scattering 
function curve of a single aggregate, which saturates for q 
values smaller than R^^, in the so-called Guinier regime 
|]l5| , is not valid down to q values of order £^^, where 
inter-aggregate correlations start to has some influence. 
Recently, a hard sphere model have been introduced by 
Posselt et al. to describe the packing of the connected 
clusters. However, their model does not take account 
of the fact that neglecting other scattering contributions 
such as scattering by thermally activated fluctuations, 
should cause S{q) and I{q) to vanish for q < £~^, since, 
for distances larger than £, the system becomes homo- 
geneous and should no longer scatter the incident beam. 
Such behavior cannot be avoided in the case of aerogels 
where £ and R are of the same order of magnitude and 
should be replaced by ^. If the experiments can be per- 
formed down to sufficiently small g-values, if no other 
artifactual inhomogeneities are present and if scattering 
by thermal fluctuations is small enough, all the experi- 
mental I{q) curves for gels or aerogels should exhibit a 
maximum. In practice, a maximum is observed, or not, 
depending on the range of g- values available. When it ex- 
ists, this maximum is more or less pronounced, depending 
on preparation, i.e. catalysis conditions. 

This paper is a review of a series of recent numeri- 
cal studies which has been, in part, published elsewhere 
0-^. We first study the larg e-q regime which is re- 
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lated to short interparticle distances within the aggre- 
gate. Then, we present quantitative results for the loca- 
tion and shape of the maximum of the I{q) curve at small- 
q values. Finally, the effect of particle size polydispersity 
in the cross-over between the fractal and Porod regime is 
discussed. The scattering functions S{q) and I{q) were 
obtained from numerical simulations of the full aerogel 
structure which is modelled by diffusion-limited cluster- 
cluster aggregation model (DLCA) in a box |^,^,^ with 
a sufficiently large initial concentration, rather than the 
hierarchical model ||2^ which is only able to build a single 
aggregate. 

2. Constraints on Theory 

2.1. The Model 

We have considered a three dimensional off-lattice ex- 
tension of the original cluster-cluster aggregation model 
po| , pT| in the case of a sufficiently large particle concen- 
tration to get a gelling network at the end of the pro- 
cess. Such a model was previously proposed by Kolb and 
Herrmann Q to describe the formation of gels, but, at 
this time, they considered a two-dimensional model on 
a lattice. Initially, identical spherical particles of unit 
diameter are randomly disposed in a cubic box of edge 
length L (here L is not necessary an integer) using a stan- 
dard sequential addition procedure: attempts are made 
to center particles, one after another, at points whose 
coordinates are random numbers uniformly distributed 
between and L. If a particle overlaps a previous one, 
it is discarded and a new trial is made. If the process 
generates N particles, the dimension-less concentration, 
or volume fraction, c is given by: 



_n N 



(1) 



From previous studies [Q, it is known that, with this 
procedure, c cannot exceed an upper limit which is called 
the "jamming concentration", Cj = 0.385. Let us con- 
sider the starting configuration as a collection of aggre- 
gates containing one particle each. At a later time, one 
obtains a collection of Na aggregates, the i-th aggregate 
containing rii particles, so that: 



(2) 



The algorithm proceeds as follows. An aggregate i is 
chosen at random according to a probability, p„. which 
depends on the number of particles rii that it contains 
given by: 



(3) 



In most of our simulations, we have taken a = —0.55, 
a value close to — where D ~ 1.78 is the fractal di- 
mension of the resulting aggregates built in 3d in 
order to insure that the diffusion coefficient of the ag- 
gregates varies with the inverse of their radius. Then 



a space direction is chosen at random among the six 
directions ±1,±1,±1 and an attempt is made to move 
the cluster by a step of one unit length in that direction 
(note that this choice corresponds to perform a random 
translational brownian motion on a lattice but, since the 
original coordinates of the particles are not integers, the 
aggregates themselves are built off lattice). If the cluster 
does not collide any other cluster during this motion, the 
displacement is performed and the algorithm goes on by 
choosing again a cluster at random, etc... If instead a 
collision occurs, the cluster is translated in the chosen 
direction by the shortest distance insuring that one of 
its particles becomes tangent to one particle of the col- 
lided cluster. Then the collection of clusters is updated: 
the two colliding clusters are discarded and a new clus- 
ter, formed by sticking together the colliding clusters, is 
added to the collection. After that, one cluster is chosen 
at random, etc... Periodic boundary conditions are used 
at the edges of the box. The process is stopped when a 
single aggregate is reached. If the concentration is larger 
than a characteristic gel concentration Cg , the final aggre- 
gate spans the box from edge to edge in the three space 
directions. This is the usual convention to define a gel 
network. 

A series of calculations has been done to determine the 
gel concentration Cg as a function of the box size L. In 
practice we have varied the concentration and we have 
performed twenty independent runs for each concentra- 
tion. The gel concentration has been defined as being the 
concentration at which ten runs end up with a gel. The 
results are given in figure 1 as a log-log plot of Cg versus 
L. As already found by Kolb and Herrmann ||23| , the 
gel concentration tends to zero in the infinite L asymp- 
totic limit. Since, at the geUing threshold, an aggregate 
of fractal dimension D reaches the size L of the box, one 
should have: 
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(4) 



Our data are well fitted with the slope —1.28 ± 0.05, giv- 
ing D ~ 1.72 ± 0.05, a value shghtly smaller but quite 
close to the fractal dimension D = 1.78 of DLCA at the 
gel point in three dimensions j2^ . 

In all the calculations reported below, we have consid- 
ered c values much larger than Cg, in order to get a gel far 
above the gelling threshold and to be sure that the cor- 
relation length ^ is smaller than the box size L. Typical 
examples of gel networks for two different concentrations 
are depicted in figure 2a and 2b where we have visualized 
a two-dimensional projection of a slice of the box. The 
slice thickness has been chosen to be proportional to i in 
order to get the same mean coverage in projection. Thus 
the only difference between the two pictures is in the way 
the apparent local densities deviate from the mean. As 
one sees larger holes in case (a), it is apparent that the 
mean cluster size is larger in that case where the con- 
centration is smaller than in case (b). For comparison. 
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we show in figure 2c a micrograph of a colloidal aerogel 
sample. 

2.2. The correlation junction g(r) 

As usual the two points correlation function (?(~f*) is 
defined such as g{~r')<fir is proportional to the proba- 
bility of finding a particle center in a volume d^r at a 
distance ~r* from a given particle center. Consequently, 
for an isotropic material, the number of particle centers 
dn located between r and r + dr from a given particle 
center is proportional to g{r)A-Kr^ dr . Knowing from eq. 
1 that, in average, the number of particle centers per 
unit volume is — , one can normalize g{r) to unity when 
r tends to infinity, by writing: 



6c 

dn = — g{r)ATTr'^dr — 2Acg{r)r^dr 

TT 



(5) 



We have used this formula to compute g{r) in the gels 
resulting from our simulations. For each particle in the 
box, we have counted the number of particle centers lo- 
cated between spheres of radius r and r + 5r taking care 
of the periodic boundary conditions when investigating 
regions outside of the box. Then we have averaged the 
result over all the N particles in the box and divided it 
by 24:cr^5r. Moreover, g{r), obtained that way, has been 
averaged over a large number of independent simulations. 
2.3. The scattering functions S(q) and I(q) 
From the single scattering theory, the structure factor 
S(ci) of a macroscopic system containing identical parti- 
cles with mean volume fraction c is given by [E6|: 



S{q) = 1 



6c 



{g{r) - l)!^Wcir (6) 
qr 



The infinite boundary of the integral means that the in- 
cident beam is scattered by the particles located within 
the whole macroscopic volume. The presence of g{r) — 1 
means that one has subtracted the intensity scattered 
by a quasi-infinite homogeneous object having the same 
boundaries as the considered macroscopic system. As a 
consequence, with formula (6), S{ci) when g ^ 
(because the q ~ contribution of the boundaries has 
been suppressed) . Quantitatively, this g = limit results 
from the following sum rules: 



6c /"^ 

— / g{r)AiTr'^dr = M - 1 
71" Jo 



and. 
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Anr'^dr = M 



(7a) 



(7b) 



where M is the total number of particles contained in the 
macroscopic volume. 

When applying formula (6) to real systems, it should 
be remembered that the underlying theory considers 
only the single scattering by static (quenched) particles. 
Therefore the scattering by thermally activated (and cor- 
related) motions of these particles is not considered. It is 



known that in the case of liquids, the latter contribution 
(which there corresponds to scattering by thermally acti- 
vated density fluctuations) gives a non zero contribution 
ai q — Q proportional to the product of temperature, bulk 
density and isothermal compressibility [ p7| . Since in the 
following, the formula will be applied to solid aerogels, 
of quite low compressibility, we will assume that such 
contribution is negligible. 

In practice, we have numerically calculated S{q) from 
the preceding g{r) results, by replacing formula (6) by: 



S{q) = l + - 



' / / X X smgr . 9 , 
{9{r)-gQ)—^^T:r^dr (8) 
qr 



Here, go is a parameter which is very close, but not 
strictly equal, to one, is an upper cut-off and the inte- 
gral is numerically computed as a discrete sum. We have 
chosen = ■§ , to avoid the boundary artifact mentioned 
above but also we have averaged g{r) over many simu- 
lations and we have limited ourselves to concentrations 
sufficiently greater than the gel concentration to obtain 
a significant range of r values below rm where g{r) ~ 1. 
Nevertheless, the truncation of the integral implies that 
the numerical results for S{q) become meaningless for q 
values smaller than: 



(9) 



and, indeed, we have observed that for q < q-^i-^ the pre- 
cise shape of the computed S{q) curve depends on both 
rm and go. To avoid this problem, we have forced the 
sum rules (7a) and (7b) to be verified, and, instead of 
using go = 1, we have computed this parameter from: 
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JJ^"' iTTr-^dr 



(10) 



where, to compute the integrals, we have used the same 
discretisation as in formula (8). We have checked that 
go is always equal to 1 within less that 0.001, however 
using (10) instead of go = 1, insures that S{q) 0, ex- 
actly, when q — > 0. This trick, which allows to obtain a 
continuation of S{q) below g^j^jj^, is expected to give a 
reasonably correct result if, in the corresponding infinite 
system, g(r) is supposed to stay constant for r > r^. 
This has been checked a posteriori by verifying that the 
numerical results are the same (within the numerical un- 
certainties) for different L- values. 

For systems formed by monodisperse particles, the 
scattered intensity /(g) can be obtained by using: 



/(g) = S{q)P{q) 



(11) 



where P{q) is the normalized form factor for spherical 
particles of unit diameter: 



P(g) = (24 



sm f - f cos f 



(12) 
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If the system contains polydisperse particles one can no 
longer calculate the scattered intensity I{q) as a product 
of P{q) and S{q). One should go back to the calculation 
of the scattered amplitude pa, which is proportional to: 



(13) 



where fi refers to the center of the i-th particle and x 
refers to a running point inside the volume of the i-th 
particle with respect to its center. The integral inside the 
sum, which should be performed over the volume of the 
i-th particle, can be calculated as a function of a^, assum- 
ing isotropy and homogeneity inside the sphere, leading 
to: 



with: 



sin(2^)-(2fi)cos(^) 



(14a) 



(14b) 



Then, assuming a random orientation of the aggregate 
over the direction of q, the scattered intensity I{q) = |v4p 
can be written as: 



/(g) = ^ A, A 



sin(gr,j-) 
qrij 



(15) 



where = {ri — fj\. 

Note that, since the i and j dependent product A,Aj 
appears inside the sum, the result cannot be split in two 
parts. In particular one cannot use the correlation func- 
tion g(r) to calculate an intermediate structure factor 
S{q). Here the double sum should be calculated directly. 
We have used such formula to calculate I{q) for single 
aggregates built using a hierarchical model. 

3. Results 

In figure 3 we show typical g{r) curve that result 
from an average over 20 simulations with L = 57.7 and 
c = 0.05. In this figure one observes a strong peak at 
r = 1, a discontinuity at r = 2, and a weaker singularity 
at r = 3 (discontinuity of the derivative). At large r, g(r) 
goes through a minimum and becomes very close to one 
(in order to avoid artifactual correlations due to the pe- 
riodic boundary conditions, we have used c values bigger 
than Cg to insure that the the homgeneous regime (when 
g{r)=l) can be reached for r < -j). For small c val- 
ues, one can observe the fractal regime, at intermediate 
r values, where g{r) follows the power law: 



(16) 



The range of this fractal regime decreases as c increases. 
Note that in the inset in figure 3, where c = 0.1, it is 
almost inexistent. 



The strong peak at r = 1 is due to the non-zero pro- 
portion of distances r = I corresponding to bonds be- 
tween contacting particles. The non-zero value of g{r) 
for r = 1+ and the discontinuity at r = 2 can be under- 
stood if one considers that g{r) can be written as: 



lir) = gi{r) +92{r) 



(17) 



where gi (r) is the contribution of couples of particles that 
are tangent to the same third one and where g2{'i') con- 
tains all the other contributions. We have observed that, 
while g2{T) is continuously varying from r = 1, where 
52(1) = 0, up to the largest distance, going through a 
maximum around r = 2^ gi{r) exists only between r = 1 
and r — 2 and reaches non-zero values at both limits. 
The weaker singularity at r = 3 can be explained by an- 
alyzing distances between spheres tangent to each sphere 
of a dimer in a manner very similar to the discontinuity 
at r = 2. 

In the inset in figure 3 we show the g{r) curves obtained 
for L — 57.7 and different c values by emphasizing the 
region near the minimum. One can see that the loca- 
tion of the minimum strongly depends on c. The mini- 
mum corresponds to distances between particles located 
at the periphery of the clusters where the local density is 
smaller. Therefore the value of r at the minimum gives 
a good estimate of the mean cluster size, let us call it f . 

The dependence of ^ with c is reported in the log-log 
plot of figure 4. Assuming that the clusters are fractal 
with fractal dimension _D, one should get: 



(18) 



The value ^ = 3 reported for the largest concentration 
c = 0.1 should be considered as artifactual since for such 
concentration the minimum of g{r) sticks on the singular- 
ity at r = 3, as seen in the inset of figure 3. A straight line 
fit of the remaining points for the largest box size gives a 
slope of -0.77 ± 0.03, leading to D = 1.70 ± 0.06, again 
slightly smaller but quite close to the fractal dimension 
D ~ 1.78 of 3d DLCA aggregates. We point out that 
in this model D increases when c increases as showed 
elsewhere [^,^. However, the above value D can be 
understood by analyzing the time evolution during the 
aggregation process at the instant where the system can 
be considered as a gelling network the fractal dimension 
of the clusters is D « 1.78 and the characteristic length 
^ value does not change compared with to measured ^ 
value at the end of the aggregation process . 

Typical S{q) curves are reported in figure 5 for differ- 
ent concentrations and for L — 57.7. The location of 
9min indicated by the arrow. All the curves exhibit 
the same damped oscillations at large q. In this figure, 
one observes a large minimum at about g ~ 4 followed by 
damped oscillations. The oscillations can be attributed 
to the 5-peak of g{r). If one considers only the distance 
contributions of the delta peak the correlation function 
g{r) writes: 



5iW = 24^'^(^ " 1) 



(19a) 
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and its Fourier transform: 



'5'oo('?) = 1 



smq 

q 



(19b) 



where z is the coordination number of the system. In 
DLCA aggregates z — 2 and the corresponding Soo{q) 
curve is depicted by the dotted hne in figure 5. As ex- 
pected, this approximation corresponds to the asymp- 
totic large-g hmit of S{q). However the large minimum 
at q ~ 4 is not accounted for by this contribution only. 
Not only the delta-peak at r = 1 influences the shape of 
the first minimum oi S{q), but also the others short range 
features of the system, as the discontinuity at r = 2 and 
the singularity at r = 3, play some role. 

The linear fractal regime is quite narrow and is more 
extended for low concentrations. It corresponds to a frac- 
tal dimension D 2± 1.7. We point out that in a system 
of homogeneous connected fractal clusters of finite size, 
the fractal dimension is only an "apparent fractal dimen- 
sion" as showed elsewhere |2^j3^. At lower q values, 
all the curves exhibit a maximum. The location of the 
maximum, qm, as well as the intensity of the maximum, 
S{qm), have been reported as a function of c in table I 
together with the corresponding values of ^. If one for- 
get the artifactual situation c = 0.1, the product q„j^ is 
almost constant, so that one has approximately: 



2.75 



(20) 



On the other hand one observes that the intensity of the 
maximum S{qm) is roughly proportional to c^^, which is 
the number of particles contained in a sphere of diameter 



S{qn 



(21) 



For very low g- values (smaller than q^) we get a linear 
behavior with slope +2. This is consistent with the low-q 
expansion of formula (7), which predicts: 



Siq) 



with: 
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C - — igo-g{r))ATrr^dr 



(22a) 



(22b) 



This result is derived from a Taylor expansion of sin qr/qr 
inside the integral of (8). Note that, for an infinite system 
(where r„j — oo and go = 1), such procedure is mathe- 
matically justified if g{r) — 1 tends to zero more quickly 
than any power law when r tends to infinity. Our nu- 
merical results, which leads to a quite size-independent 
result for C, strongly support a large-r exponential decay 
for g{r) — 1. This defines another characteristic length ( 
which can be usefully compared to ^. In table I, we have 
reported the numerical estimates of ^ from our computed 
S{q) curves for the different c values considered. We find 



that C is roughly proportional to ^, and therefore ( does 
not bring any new information on the structure. 

As it should be expected, the characteristics of the 
maximum of S{q) do not depend only on the mean clus- 
ter size f , but also on the extension of the cluster size 
distribution. The cluster size distribution is reminiscent 
of the size distribution of the aggregates observed dur- 
ing the aggregation process and therefore it strongly de- 
pends on the nature of the aggregation mechanism. As 
a consequence, all the quantitative analysis done above, 
as well as the values of the constants appearing in for- 
mulae (20) and (21), are only valid for DLCA where it 
is known that the aggregate size distribution presents a 
well defined maximum. 

When the cluster perform a ballistic random mo- 
tion instead of diffusion brownian during the aggre- 
gation process described in section 2.1., one obtains 
the ballistically-limited cluster-cluster aggregation ^] 
(BLCA) model which is able to describe the aerosols 
in the Knudsen regime. In BLCA the polydispersity of 
cluster sizes is larger than in the DLCA model. On the 
other hand, when a very small sticking probability is in- 
troduced in the DLCA model one obtains the chemically- 
limited cluster-cluster aggregation |3^j3^ (CLCA) model 
where it is known that the aggregate size distribution is 
broader than the DLCA and BLCA models. In figure 6 
we show that this polydispersity features are observed in 
the S{q) curves, where the broadest peak at q^ corre- 
sponds to the CLCA model. 

The log-log plots of I{q) versus q for different concen- 
trations are reported in figure 7. As a result of the mul- 
tiplication by P{q) and the use of the logarithmic scale, 
the maximum appears to be relatively less pronounced 
than in the S{q) curves. The large-q oscillations in the 
Porod region are due to the monodisperse particles form- 
ing the system. These oscillations become more and more 
damped when one increases the particles polydispersity 
jlSt as it will be showed in figures 10. 

4. Discussion 

In this section we would like to discuss the theoretical 
results in the light of experimental data. 

4-.1. Short-range correlations (large q-values) 
Colloidal silica aerogels have been prepared using the 
process described in ref. [10[ . They have densities ranging 
from 0.070 to 0.380 gjcnr. They are made of small col- 
loidal spherical particles with a quite low diameter poly- 
dispersity as it has been checked on electron micrographs 
such as figure 2c. We have checked that the intensity 
scattered by the correpo ndly d iluted sol can be fitted by 
an averaged form factor P{Q) given by: 



PiQ) 



P{Q, a)g{a)da 



(23a) 



where P{Q^ a) is the form factor of a spherical particle 
of diameter a and g{a) a truncated gaussian distribution 
given by: 



(23b) 
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The same analysis has been done for all our experimen- 
tal data and in the following we give the res ults for the 
structure factor as being the ratio of I{Q) by P{Q) as a 
function of the reduced wavevector q — Qa^. All the ex- 
perimental S{q) curves, determined this way, have been 
normalized such that S{q) — > 1 for q —>■ oo. It has been 
shown [ p^t that using formula (11) in presence of small 
polydispersity, in spite of this approach, the S{q) curves 
present the same broad minimum. 

In figure 8a one gives log S{q) versus logg for different 
aerogels made of particles of the same size (ao = 96A ) 
but with densities ranging from 0.070 to 0.250 g/cm?. 
On this figure, one observes the same characteristic broad 
minimum followed by damped oscillations that we have 
observed in the simulated curves, and the density fixes 
the size ^ of the clusters. In figure 8b we compare two 
experimental S{q) curves for the same aerogel density 
{p = O.lOg/cTO^) with the simulated curve in the diffu- 
sion limited case. It is clear in this figure that the agree- 
ment between theory and experiments is only qualitative. 
Even if the data are very noisy for large q values it seems 
that the large q oscillations of the experimental curves 
are more damped. But we would like to focus on the 
larger discrepancy which is that the minimum is wider 
and deeper in the experimental curves. This discrepancy 
is systematically more important for bigger particles, and 
its cannot be attributed to the kind of restructuring as 
showed in ref. |l7| because a restructuring mechanism 
in DLCA aggregates increases its coordinations number 
resulting in the S{q) curves a deeper minimum but not 
wider. Since it is difficult to imagine some other realistic 
restructuring mechanisms able to fully account for the 
observed discrepancies, we do not trust the earlier inter- 
pretations which were considering quite large coordina- 
tion numbers One might invoke other possible 
explanations for the discrepancies such as small-g modifi- 
cations of the form factor or corrections to the scattered 
intensity due to some shape deformation of the particles 
near their contact zone. This last effect might be ap- 
proximately taken into account by considering a differ- 
ent length for the particle-diameter and for the center- 
to-center distance between contacting particles. However 
all these considerations, if they might sometimes give a 
better fit, appear to be too "ad hoc" to really improve 
the comprehension of the problem. 

Here we would like to propose another tentative inter- 
pretation. In general a complete theory of scattering (in- 
cluding multiple scattering, shadowing, refraction etc..) 
should consider two dimensionless parameters, Qa and 
ka = Sttj. The fact that the theoretical S{q) curve con- 
sidered above does not depend on the extra parameter ka 
comes from all the considered approximations. However 
the simple scattering theory should be recovered in the 
limit ka 0. Some corrections might appear for large 
ka values. The fact that in figure 8b the theoretical curve 
can be considered as the limit of the experimental ones 
when a — > support this analysis. Moreover the param- 
eter ka is quite large in our case. We used a combination 



of two incident neutron wavelengths of 6A and 18A in 
the experimental set-up. Thus our ka values are in the 
range 30 to 300, close to the values involved in the geo- 
metrical optics approximation. It is reasonable to admit 
that corrections to the simple scattering theory, such as 
shadowing, refraction, multiple scattering effects, cannot 
be neglected for such large values. 

4-.2. Long-range correlations (small q-values) 
As is known, colloidal aerogels are formed by particles 
bigger than standard aerogels. This fact impedes to ob- 
tain information about its long-range correlations (small 
q- values). For study the long range correlations we have 
considered SANS experiments on standard silica aerogels 
[^,|^. They are prepared by chemical reactions (hydrol- 
ysis and condensation) of organosilicates. According to 
the ph value of the hydrolysis aqueous solution, we can 
distinguish "basic" and "neutral" aerogels. The basic 
aerogels 34 1 are made of larger sized, but strongly poly- 
disperse, primary particles while, for the neutral aerogels 
particle sizes are smaller and extend down to the 
atomic scale. According to previous studies |^,^,|lO|, 
only colloidal and basic aerogels can be considered as 
grown according to DLCA while the neutral aerogels are 
more likely grown according to the chemically-limited 
cluster-cluster aggregation process 32,3^]. Three ex- 
perimental I{q) curves for basic aerogels with different 
densities are compared with simulations in figure 9. The 
concentrations that we have used for the fit are those 
corresponding to the aerogel density p, according to the 
formula: 



Po 



(24) 



where po = 2, 2g.cm~'^ is the density of silica. The only 
two adjustable parameters are the mean particle diam- 
eter value and a multiplicative constant for the inten- 
sity. Note that varying these parameters in a log-log plot 
does not change the shape of the curve but only leads 
to translations along the x and y directions. The three 
fits have been performed together using the same values 
for these parameters. A discrepancy occurs in the Porod 
region where the experimental curves stay slightly below 
the maxima of the large-g oscillations of the theoretical 
curves. Such discrepancy can be attributed to the strong 
polydispersity of the primary particles jist and will be 
discussed in section 4.3.. 

We point out that, our DLCA model neglects rota- 
tional diffusion, aggregates deformations as well as all 
kinds of restructuring effects. However, as mentioned 
above, large restructuring is certainly not present, at 
least in the early stages of the aggregation process. But, 
in the last stages, when the gel structure is under for- 
mation, the diffusion process looses its leading character 
and it might be that restructuring and rotational mo- 
tions have some influence on the gel structure and the 
inter-clusters correlations. 

4-3. Effect of the particle size polydispersity 
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To study the influence of the particle size polydisper- 
sity we have calculated I{q) , using formula (15), for hier- 
archical aggregates made of spheres whose diameter are 
distributed according to a truncated gaussian distribu- 
tion characterized by Ce/ / = ^ i where oq and a are the 
mean value and the standard deviation of their diameters 

It has been shown ||18[ that the crossover wave-vector 
Qc between the fractal and Porod regime, for (Je// > 0.1, 
varies as: 



l.GcTe//) 



(25) 



where (Jeff is the dimensionless standard deviation, or 
polydispersity. Note that, when (Jeff = the crossover 
wave- vector qc is equal to as can be verified in fig- 
ure 7. In figure 9 the experimental curves exhibits the 
crossover take place at Qc approximately equal to ^ and 
according the above equation Jeff should be equal to 
0.25. We have computed, for simplicity, a single aggre- 
gate with a particle number N equal to 32 using the 
hierarchical procedure with a particle polydispersity 
aeff = 0.25, the theoretical I{q) curve has been calcu- 
lated using equation (15) and has been averaged over 32 
simulations. The resulting curve is shown in figure 10a 
(solid line) in a log- log plot of I{q)q'^, emphasizing the 
large-g region, and for comparison we have depicted the 
case when aeff — (dashed line) and the experimental 
curve for c = 0.043(open circles). We can note in this fig- 
ure that introducing a polydispersity not only the large-g 
damped oscillations disappear, but also qc shifts towards 
low-g values. This result can be attributed to the fact 
that larger particles dominate the scattering. Even if a 
single aggregate does not account for the realistic long- 
range correlations, figure 10a shows the good agreement 
at large g-values between the theoretical I{q) curve for 
Jeff = 0.25 and the experimental ones. Figure 10b shows 
a log-log plot of I{q) versus q. Here the dashed line rep- 
resents the numerical I(q) results at the low-(g) values 
in a system of connected fractal clusters for c = 0.043, 
and solid line represents results obtained at large-g values 
of a single aggregate made of polydisperse particles. In 
this figure the agreement between the simulation and the 
entire experimental results (open circles) is remarkable. 
Conclusion 

In this paper we have shown that the diffusion-limited 
cluster-cluster aggregation model (DLCA) can explain 
the structure of "colloidal" as well as "basic" aerogels, as 
revealed by small-angle neutron scattering experiments 
since in most cases the full I{q) curve can be quan- 
titatively accounted for. This is a real progress com- 
pared to previous approaches which were focusing on the 
intermediate-g fractal regime only. Such modelization 
is now used to study some other physical properties of 
aerogels. In the present issue we give two examples of 
application which are numerical calculations of aerogel 
sintering 1 35 1 and simulations of phase transitions in the 
pores of aerogels pq] . We also present a numerical study 



of the evolution of the I{q) curve during the aggregation 
process Other applications are under progress. Fur- 
ther investigations are needed to account for the experi- 
ments on "neutral" aerogels which exhibit a larger fractal 
dimension than DLCA aggregates. Neutral aerogels are 
made of very small and strongly polydisperse particles 
whose sizes extend down to the atomic scale. It might 
be that their growing mechanism is closer to chemically- 
limited aggregation (CLCA) than DLCA but also there 
might exist some complex restructuring mechanisms due 
to their flexibility. 

One of us (A. H.) would like to acknowledge support 
from CONICIT (Venezuela). 
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C i qrn gm| S{qm) C £ 

0.025 10.7 0.26 2.78 24.7 0.81 38.5 3.6 

0.033 8.4 0.31 2.60 16.0 0.82 28.1 3.3 

0.038 7.8 0.37 2.88 12.8 0.71 25.9 3.3 

0.043 6.7 0.40 2.68 10.5 0.81 21.0 3.1 

0.050 6.3 0.43 2.71 8.3 0.66 18.4 2.9 

0.061 5.3 0.52 2.76 6.2 0.68 15.5 2.9 

0.100 3.0 0.80 2.40 3.3 1.22 10.0 3.3 



TABLE I. For each concentration c considered in the simulations, we have reported the location ^ of the minimum of g(r), 
the location qm of the maximum of S{q), the product qm£,, the intensity of the maximum S{qm) and the ratio ^^^y-* , the 

characteristic length entering the low-q expansion of S{q) and the ratio |. We recall that q is here a dimensionless quantity 
which is, in fact, equal to 2Qro where Q is the dimensioned wavevector and ro is the radius of the primary particles and also 
that S{q) has been normalized to unity for large q. 
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FIG. 1. Log-log plot of the gel concentration Cg as a function of the box size L. The fit by a straight line is shown which 
gives a slope of —1.28 ± 0.05. 

FIG. 2. Two dimensional projections of the particles contained in a slice of thickness i after obtaining a gel in a box of size 
L = 57.7. Cases (a) and (b) corresponds to c = 0.0095, i = 34.6 and c = 0.038, i = 8.65, respectively. Cross sections of particles 
that are cut by the front slice edge are shown in black, (c) Micrograph of a 270 A colloidal aerogel sample. 

FIG. 3. Plot of g{r) versus r for L = 57.7 and c = 0.05. Inset: g(r) curves for different c values (c = 0.025, 0.05, 0.1). These 
curves result from averages over 20 simulations. 

FIG. 4. Log-log plot of ^, location of the minimum of g{r), versus c, for L = 28.8 and L = 57.7. The fit by a straight line of 
the L = 57.7 data, excluding the c = 0.1 point, is shown giving a slope of —0.77 ± 0.03. 

FIG. 5. S{q) versus g for L = 57.7 and different c values: c=0.025, c=0.05, c=0.1. The dotted line correspond to Soo{q)- 
Here and in the following figures the location of 9^;^ is indicated by the arrow. 



FIG. 6. Log-log plot of S{q) versus q for L = 57.7 and c = 0.05, the solid line, dashed line and dotted line, correspond to 
DLCA, BLCA and CLCA model, respectively. These curves result from average over 20 simulations. 



FIG. 7. Log-log plot of the scattering intensity curve I{q) versus q for different c-values. The parameters are the same as in 
figure 5. 

FIG. 8. (a) Experimental S{q) curves (obtained as explained in text) for the 96A -colloidal aerogel family. Samples are 
labeled by their densities, (b) Experimental S{q) curves for samples of various particle diameters but with the same density, 
p = O.Wg/cm^, and the simulated S{q) curve in the diffusion-limited case. 

FIG. 9. Comparison between simulations and experiments for three base catalyzed aerogels of the same family. The concen- 
trations used in the simulations c = 0.033, 0.043, 0.05 are calculated from the aerogel densities (p — 0.073, 0.095, O.llOg.cm"^) 

using forunila (11). The two adjustable parameters, which arc a multiplicative constant for the intensity and the particles 
diameter, taken to be 46A, are the same for the three curves. The curves have been arbitrarily shifted vertically for clarity. 

FIG. 10. (a)Log-log plot of I{q)q* vs. q for (Te// = and c = 0.043(dashed line), and for iTe// = 0.25 with = 32 (solid 
line). The open circles represent the experimental curve for c = 0.043. (b) Same data as in (a) but with a more extended q 
interval and represented in a log-log plot of I{q) vs. q. 
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Fig 10b 




